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Optical imaging methods are typically restricted to a resolution of order of the probing light 
wavelength \ p by the Rayleigh diffraction limit. This limit can be circumvented by making use of 
multiphoton detection of correlated A^-photon states, having an effective wavelength \ P /N. But the 
required A-photon detection usually renders these schemes impractical. To overcome this limitation, 
recently, so-called optical centroid measurements (OCM) have been proposed which replace the 
multi-photon detectors by an array of single-photon detectors. Complementary to the existing 
approximate analytical results, we explore the approach using numerical experiments by sampling 
and analyzing detection events from the initial state wave function. This allows us to quantitatively 
study the approach also beyond the constraints set by the approximate analytical treatment, to 
compare different detection strategies, and to analyze other classes of input states. 
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I. INTRODUCTION 

Conventional optical lithography schemes are impeded 
by the Rayleigh diffraction limit P, Q owing to the wave 
nature of light. Not least due to the continuous and 
increasing demand of faster devices in microelectronics, 
smaller structures in semiconductors need to be written. 
But highly energetic photons can lead to enhanced ma- 
terial damage and demand the development of new op- 
tical imaging systems. Motivated by this, several alter- 
native routes have been explored. A prominent example 
is quantum lithography [3|-|6| based on number-path en- 
tangled states of light. But generation of the required 
entangled states, decoherence due to photon loss and de- 
tection remain a challenge. Alternatively, classical light 
fields can be used. One tool is post-selection of certain 
detection events which contain image data with high res- 
olution 0, HJ- Other methods aim at facilitating a TV- 
photon absorption process leading to a reduced effective 
wave length [9l4l3l|. These, however, require multiphoton 
transitions, which complicate a practical implementation. 
This requirement can be relaxed by means of dark states 
in multi-level svstems [3] or of Rabi oscillations in two- 



level systems [15 1 



Recently, a different method has been proposed to re- 
move the requirement of multi-photon detection [l|| . In 
the so-called optical centroid method (OCM), the N- 
photon detector is replaced by an array of single-photon 
detectors. For each input pulse, the results of the indi- 
vidual single-photon detectors are then combined into a 
single variable, the centroid. The probability distribution 
of the centroid obtained from many repetitions of the 
measurement then contains the desired sub-wavelength 
resolution. Thus, no A-photon absorber is required, and 
at the same time, the detection becomes more efficient as 
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all events rather than only A-photon events contribute 
to the signal. Recently, a proof-of-principle experiment 
of OCM for two photon entangled state has been carried 
out by [It}. The experiment used two kinds of detec- 
tors: One is a pair of single photon detectors for OCM 
while the other is two-photon absorption detector for QL. 
OCM achieved the same quantum supcrresolution as QL 
but much more efficiently. 

Motivated by this, here, we study the OCM using nu- 
merical experiments. In contrast to existing approxi- 
mate analytical approaches, the numerical study allows 
us to quantitatively study OCM over a broader parame- 
ter range, to compare different detection strategies, and 
to analyze other classes of input states. In particular, 
in addition to the NOON and jointly Gaussian input 
states proposed in [ill, we analyze cat states made out 
of superpositions of coherent states. For our analysis, 
we first generate random detection events from position 
space wave function. Afterwards, we model the detection 
by a suitable discretization using different methods, and 
apply the centroid method. We characterize the error 
in the recovery of the theoretical centroid distribution 
by root mean square (rms) deviation from the original 
wave function. We confirm our method by recovering key 
predictions from the analytical calculations [lj| quanti- 
tatively, such as the dependence of multiphoton count 
rates for jointly Gaussian states on the properties of the 
input state. 

We then study the dependence of the rms deviation 
on the size of the individual detectors in the array. We 
find several features which may appear counter-intuitive 
at first sight, but can be explained in terms of the mea- 
surement procedure chosen and of statistical effects most 
pronounced in the small detector limit. For intermediate 
detector sizes, the rms deviation is nearly a linear func- 
tion of the detector size. We analyze this dependence for 
NOON-states with different A and show how the slope 
of the linear dependence can be traced back to the struc- 
ture of the underlying wave function. We then extract 
single-and multi-photon detection events, and verify the 
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improved efficiency expected for the OCM over iV-photon 
absorption schemes. For the jointly Gaussian state, we 
augment this analysis by studying the rms deviation for 
different N as a function of the detector size while keep- 
ing the feature size constant. We finally turn to coherent 
cat states, and show that the OCM can also be applied 
to them. These states are of interest, as the feature size 
of the cat states can continuously be tuned via the co- 
herent state parameter a. Furthermore, the structure of 
the cat state wave function along the centroid axis can be 
modified via the phase of a such that it differs from the 
regular patterns observed for NOON or Gaussian states. 

The paper is organized as follows. In Scc.|lH we provide 
the theoretical background. In particular, we review the 
OCM method in III A[ introduce the various non-classical 
input states in Sec. Ill B| discuss the implementation of 
the numerical experiments (Sec. Ill C[) and the detection 
system (Sec. Ill Dp . and describe the procedure to esti- 
mate the error in the OCM recovery in Sec. Ill El In 
Sec. IIIII our results for different types of non-classical 
input states are discussed. The final Sec. II VI summarizes 
the results. 



II. THEORETICAL CONSIDERATIONS 
A. Centroid method 

We briefly review the formalism of resolution enhance- 
ment developed in fill [l8j using quantum lithography 
and optical centroid measurements. We consider pho- 
tons of wavelength X p incident on a detection plane un- 
der an angle such that the wavelength associated to the 
wave vector component in the detection plane is given 
by A = Xp/ smO. In the following, we restrict the analy- 
sis to one spatial dimension in the detection plane. Let 
a(k)(a' (k)) denote the photon annihilation (creation) op- 
erator in this transverse momentum space. 

These operators follow the commutation relation 



[a(k),tf(k'j\ =6(k-k') 



(1) 



An A-photon momentum eigenvector can be defined us- 
ing a' (ki) as 



|fci,fc 2 



tf(k 1 )tf(k 2 )...a?(k N )\0), (2) 



where |0) indicates the state without photons. The mo- 
mentum space wave function representation of a pure A^- 
photon Fock state \N) is given by 



<p(ki,k 2 ,...,k N ) = (ki,k 2 ,...,k N \N). 



(3) 



The transverse momenta of thephotons are restricted by 
the Rayleigh diffraction limit [lj such that 



(j)(k 1 ,k 2 ,...,k N ) = 



(4) 



for any \k n \ > 2tt/\. The spatial annihilation operator is 
defined as 



A(x) 



1 



/2tt 



dka(k) e 



ik.t 



(5) 



A corresponding A^-photon state can be constructed as 
1 



\Xl,X 2 , 



,XNl 



A^ Xl )...A^x N )\0) 



(6) 



The position space wave function is obtained by taking 
the A^-dimensional Fourier transform of </>(fci, k 2 , . . . , k^) 
and is given by 



1p(xi,X 2 , . . .,X N ) = (xi,X 2 , ■ ■ .,X N \N). 



(7) 



Of fundamental importance in the centroid scheme is the 
transformation of the photon position coordinates xi to 
the centroid coordinate X and the relative position co- 
ordinates £„ which are defined as 



X 



1 



N 

11=1 



x. 



(8) 



Note that when x± = x 2 = ■ ■ ■ = ijy = x, X = x, that 
is, if all photons are constrained to arrive at the same 
point x, the relative positions become zero and the 
centroid variable also lies at the same point x according 
to Eq. ©. 

Conventional quantum lithography proposals rely on 
A^-photon detection at a single detector at position x. 
The input state is a A^-photon state pm = \N)(N\. The 
probability density for multi photon absorption is there- 
fore given by 

p MP (x) = {:I N (x) :) = {[A\x)] N [A(x)} N ) 

= N\ij(x,x,...,x)\ 2 . (9) 

The centroid method, in contrast, relies on an array 
of detectors, such that each of the incident A^ photons 
is detected by an individual detector. Out of the indi- 
vidual photon detection positions, the centroid position 
X is calculated. The probability density of the centroid 
distribution is obtained from 

(: < =1 /> n ) :) (10) 

as a marginal probability density by tracing out all rela- 
tive positions £„ using x n — X + , 

Pocm(X)= f dS 1 ,...,dt If - 1 (:I% =1 I(X + S n ):) 



.,d^ N - 1 N\ 1 p(X + ^ 1 



,x + a N )\ 2 . (ii) 



The key result of Tsang's analysis is that the two dis- 
tributions Eq. (jlip and Eq. (|9]) have the same spatial 
dependence and thus resolution. Therefore, the techni- 
cally challenging A^-photon detection in a single detector 
can be replaced by A^ single-photon detections. A further 
advantage of the centroid method is that it is potentially 
much more efficient than multi-photon detection, i.e., for 
experimentally relevant parameters, Pocm(X) can ex- 
ceed pmp{ X) significantly [l6l - [l8j . 
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B. Trial states 

Throughout our analysis, we apply the optical cen- 
troid method using numerical experiments to various 
non-classical photon states, which we introduce in the 
following. 



1. NOON states 

Following [lH, we consider momentum-correlated TV- 
photon NOON states 



ITVOOTV) 



1 



V2 



(\N)a\0)b + \0)a\N) b ) 



112) 



in which the two modes A and B correspond to two wave 
vector directions for the photons incident on the detec- 
tion plane, with mean transverse component of the wave 
vector ko and — ko, respectively. The momentum spread 
of the transverse wave vector is denoted by Afc with 
fco 3> Afc. Denoting the (normalized and even) trans- 
verse wave vector profile as /(<?), the momentum space 
wave function is 



'Jnoon 



(fci, fc 2 , . . . , k N ) = 



1 



V2Afc^ 

JV 

n/ 



N 



.71=1 



k n fco 
Afc 



fco 



Afc 



(13) 



By Fourier 

--l/i ID - q 2 /(2Ak 2 ) 



transform, assuming f(q) = 
, we then obtain the position space 



V2Ak N 



wave function 

iI>noon(xi,x 2 , ■ ■ ■ ,x N ) 

if ' - 

X COS I fco X{ 

\ 7=1 / 

The probability distribution thus follows as 



-iAfc 2 Ej^xJ 



TV 



(14) 



\iPnoon(xi,x 2 , . 

N 

2 1 — 1 e 



■ ,xn) 



Afc 2 Ef =1 A 



JV 



k j2- 



(15) 



Note that this expression is only normalized in the limit 
Afc <C fco, in which the field operators for the two modes 
A and B satisfy the usual bosonic commutation rules. 

For numerical simulations we employ dimcnsionless po- 
sition variables x, = x/X such that the position coordi- 
nates are measured in the units of A. In terms of the 
dimensionless variables, Eq. (|15[) becomes 



'Afc 



N 



N 



COS 



(16) 
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FIG. 1. (Color online) Momentum (a-c) and the correspond- 
ing position (d-f) space probability distributions of a two- 
photon NOON state (top row), joint Gaussian state with 
TV = 2 and B — /3 = A -1 (middle row), and correlated cat 
state with a = i (bottom row). 



Using A = 27r/fco, we have defined the dimensionless 
width a = fco/ Afc. Throughout the numerical analy- 
sis, we have chosen a = 4v2~7r such that Afc <C fco as 
required. 

An example for the position space probability distri- 
bution of a two-photon NOON state using dimensionless 
variables is shown in Fig. (Q]). 



2-photon NOON states can be created using sponta- 
neous parametric down conversion combined with the 
Hong-Ou-Mandel effect HES]. Such NOON states have 
been employed in proof-of-principle experiments on res- 
olution enhancement [f| [13, H9|- Higher order NOON 
states with larger TV have also been created experimen- 
tally mum. 
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2. Jointly Gaussian states 



states of the form \2t 



The momentum space representation of a jointly Gaus- 
sian state can be written as fl^ 



4>JG(ki,k 2 , ■ 
with components of the matrix B given by 

I'll n 

B nrr}± — 



k n B nm k„ 



(17) 



1 

N2B2 

1 



N 2 B 2 N/3 



(n m) 



(18) 



C is a normalization constant. B can be interpreted 
as the width of the average momentum of the photons, 
whereas /3 characterizes the width in momenta relative 
to the average [l8| . These two parameters determine the 
variance of k n given by 



1 



(19) 



The position wave function for the jointly Gaussian state 
then follows as 

\ip JG {xi,x 2 , ■ ■ -,x N )) oc cxp(-^ 

The multi photon absorption pattern (: I (x) :) is a 
Gaussian with a root mean square width 1/(2NB). If 
B = f3/y/~N, the distribution becomes classical and the 
photons are independent. This distribution is a symmet- 
ric Gaussian with equal widths along the diagonal axes 
in the configuration space. The classical multi photon 
absorption width is given by the standard quantum limit 
Wc = l/(2-\/A r (fc^)). Increasing B from the classical 
value makes the distribution quantum mechanical. The 
maximum value of B, B = (fc 2 ,) renders the ultimate 
quantum limit for the width of the multi photon absorp- 
tion pattern W min = l/(2^^ v /(fe2)). 

The quantum mechanical elongated jointly Gaussian 
has a shorter width along the centroid direction and is 
stretched along all the other diagonals in the position 
space. 

The position space probability distribution is 

\lp.JG{xi,X2, ■ ■ ■ ,Xn)}\ 2 OC exp(-2 ^ XnB'^Xm) ■ (21) 

n,m 

Fig. ([1]) shows the position space probability distri- 
bution for a quantum mechanical jointly Gaussian state 
in terms of dimensionless position variables defined in 
Sec. (|ITBT|) . 



3. Cat states 

The third class of states we investigate are correlated 
coherent cat states. We in particular analyze two mode 



IVccc) =Af(|a)i|a) 2 + |-a)i|-a) 2 ) ! (22) 

in which \a)i is a coherent state in mode i with com- 
plex parameter a, defined by a,i\a)i = a\a)i where Oj 
is the photon annihilation operator of mode i. Note 
that other two mode cat states have also been consid- 
ered in the literature in different contexts [2Ti - 429| . These 
states can be seen as an extension of the famous sin- 
gle particle Schrodinger coherent cat states [3(3] given 
by |a+) = | a) + \ — a) [3l|. A proposal for the gen- 
eration of Schrddinger's cat state has been discussed 
in [32j . Experiments performed on the optical cat states 
include [33-35]. Af is a normalization factor which is 
given by 1/^/2(1 + exp(-4|a| 2 ). 

While Eq. flU]) resembles a two-photon NOON state, 
there is an important difference between the states. The 
bi photon NOON state consists of exactly two photons. 
Individual realizations of the same cat state, however, can 
contain different numbers of photons. The probability 
p{n\,ri2) to detect n,; photons in mode i evaluates to 

p(ni,n 2 ) = (n 1 ,n 2 \ip CC c)(tp CCC |^1 ) ^2/ 
,e- 2 H 2 



=2AP 



nil n 2 ! 



| a |2(m+n a ) [l + (_l)«i+™2] 



(23) 



It follows that the state can only contain even numbers 
of photons. For a certain range of \a\, states with less or 
equal to two photons are populated with high probability. 
This motivates the analysis of detection events, in which 
one photon is detected in each of the two modes. In the 
following, we will focus on this case. 

The position representation of a single mode coherent 
state is 13611 



1 



(x a 



y/xoVn 



-.e 2 



k(a 2 -\a\ 2 )-^(x/x -V2a) 2 



(24) 



with xq = (muJo). The position space representation 
of Eq. then follows as 



ipccc(xi,x 2 ) =(xi|a)i(a; 2 |a) 2 + (jei| - a) 1 (x 2 \ - a); 

^ -^z( x i+ x l) - 2 |a|acos(^) 



x cosh \ ^^-(xi +x 2 ) \ , 



(25) 



where a = \a\ exp(i4>). 

For the further discussion, we now specialize to the 
case of a = i\a\ (i.e., <j) = n/2) to obtain a probability 
density 



\ip ccc (xi,x 2 ,a = z|a|)| 2 
4W 2 -\(xl+xl) 



'V2\c 



nxZ 



•'•() 



-{xi + x 2 ) 



(26) 
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The structure of this expression resembles the probability 
density for the two photon NOON states. Other choices 
for <p allow to rotate the fringe pattern in the position 
space with respect to the centroid axis. This way, the 
condition for the centroid method that the probability 
density should separate into a function of the centroid 
coordinate times a function of all other coordinates can 
be continuously violated by modifying <f> from 7r/2. 

As before, for the numerical calculations we define 
dimensionless position variables x. For x = 2t:xqx 1 
Eq. (|26|) becomes 



AJ\f 2 



i a 



e -i^(xi+xi) cog 2 



(2irV2\a\(xi+X2)) . (27) 



\1>c 



Comparing Eq. ([27]) with Eq. (fT5j) . we find that the ex- 
ponential parts of the two expressions agree for a = 1, 
whereas the cosine-parts differ in the sense that the cat 
state allows to continuously tune the structure size of 
the fringe pattern by varying \a\. This way, the sub- 
wavelength resolution capabilities of the centroid method 
can be probed independent of the number of detected 
photons. 

Note that alternatively, one could choose dimensionless 
variables x = y/2n xox/\a\, such that the cosine-parts of 
the NOON-state and the cat-state probability densities 
become the same. However, in this scaling, neither the 
overall scaling nor the width a of the cat state are inde- 
pendent of the choice of \a\. 

The probabilty distribution for <fi = 7r/2 in configura- 
tion space is shown in Fig. ([T]). Increasing <j> from 7r/2 
(Fig. (jTJ) results in an increase of the width of each fringe 
along the centroid direction such that the fringes merge 
into each other. Further increase causes the middle fringe 
to become smaller and smaller until eventually at <f> w w, 
the two modes of the coherent states no more interfere 
and can be seen as separate circles. 



C. Numerical experiments 

In order to simulate the centroid method numerically, 
we first generate a large number of random photon events 
distributed and correlated according to the respective po- 
sition space wave function. Next, we apply a discretiza- 
tion scheme to model the experimentally accessible signal 
for different detector characteristics. Finally, we perform 
the centroid analysis on this data to recover a wave func- 
tion. We then compare this wave function to the origi- 
nal position space wave function and calculate the root- 
mean-square deviation to access the predictive power of 
the centroid method for the given conditions. In the fol- 
lowing, we describe all steps in more detail. 



1. Random number generation 

The first step is to generate random photon events ac- 
cording to the respective position space wave function. 
As can be seen from Eqs. ([T5|). and ([2"7). the po- 
sition variables are correlated. To facilitate the random 
event generation, we apply unitary transformations to a 
set of uncorrelated variables. For example, the probabil- 
ity distribution 



Afc 
0? 



-Afc 2 (xl+xl) 



[ko(xi + x 2 )] 



(28) 



of Eq. (|15p specialized to the case of two photons (N = 2) 
does not allow to draw random positions x\, x 2 indepen- 
dently. Applying the transformation 



Y 2 — M 2 X 2 , 

with X 2 = (x 1 ,x 2 ) T , Y 2 = (yi,y 2 ) T , and 

1 



(29) 



Mo 



1 -1 
1 1 



the probability distribution Eq. 

2 



becomes 



Ak 



-Ak'y 



00 



V2k. 



0V2 



(30) 



In this form, random values can be drawn for y\ and y 2 
independently. In a similar way, also 3- and 4-photon 
events can be handled, with coordinate transforms 



1 

71 



1 1 1 

y/3/2 -a/372 
a/2 -1/V2 -1/V2 



and 
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1 


1 








V2 


-V2 


7t 





2\/273 


a/273 


-a/273 




_V6 


-1/V3 


-1/V5 


-I/a/3. 



M 4 



respectively. 

After application of the variable transform, for the 
NOON and correlated coherent cat states, N — 1 of the 
new variables follow univariate normal distributions. The 
only exception is one variable which coincides with the 
centroid coordinate except for an overall scaling factor 
of \ N . For the random number generation in this vari- 
able, we apply the cumulative distribution function. This 
function for a continuous random variable V is obtained 
from the probability density function T via 



X(t)dt. 



(31) 



If random numbers drawn from the cumulative 

distribution, J-~ 1 (xi) is a random sample from I. 
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then follow from Eq. ([8]) as 



i 



2d 
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FIG. 2. (Color online) Model for the detection setup, (a) 
shows an array of detectors of size do- The red arrows indicate 
the central positions of the individual detectors, (b) shows the 
corresponding array with double detector size do. Compared 
to (a), only half of the possible detection positions occur, (c) 
shows the array of (b) suitably shifted that now the possible 
detection events cover those of (a) which are missing in (b). 



The situation is simpler in the case of a jointly Gaus- 
sian distribution. In this case, all transformed variables 
follow univariate Gaussian distributions, even though one 
of the new variables is proportional to the centroid vari- 
able. 

A back-transformation from the variables j/j to the 
original coordinates x% via the inverse M i 1 then yields 
the desired correlated photon events. In all cases, we ver- 
ified that the position distribution of the obtained ran- 
dom numbers agree to the original respective position 
space distribution functions. 



D. Modeling of the detection system 

The positions of the correlated iV-photon states ob- 
tained in the previous section are continuously dis- 
tributed. Any measurement, however, employs detectors 
of finite size, which yield a discretized position informa- 
tion. To model this discretization, we assume an array 
of identical detectors, as shown in Fig. [21a). Each de- 
tector has size do- We assume that the detectors do not 
overlap, and that there is no space between two adjacent 
detectors. Thus, each photon position of the continuous 
position distribution obtained from the random number 
generation can uniquely be assigned to a single detector. 
The numerical discretization procedure thus amounts to 
replacing the position of each individual photon with the 
central position of the detector which registers the re- 
spective photon. 

In a suitable coordinate system, the possible measure- 
ment outcomes are ido, with i € Z. The possible out- 
comes for the centroid coordinate in a iV-photon detec- 
tion with photons hitting the i n th detector (1 < n < N) 



X 



do 



N 

N 

n=l 



(32) 



If the detector size is changed, also the set of possible 
detection outcomes is changed. For example, in Fig.[2jb), 
detectors of double the size compared to those in (a) lead 
to a set of possible measurement outcomes which contain 
only every second possibility compared to that of (a). 
In addition, the detection array can be shifted, with an 
example for shift do shown in Fig. [5Jc). The measure- 
ments of (b) and (c) together allow for a set of potential 
measurement outcomes which coincides with that of (a). 
Note, however, that the detectors in (b) and (c) overlap 
spatially, such that the combined measurements of (b) 
and (c) are not equivalent to measurement (a). 

Throughout the later analysis, we will compare the 
performance with different detector sizes. In order to 
obtain comparable predictions for the wave function, we 
proceed as follows. For the smallest detector size do, we 
apply the centroid method for a single detection array 
position, which we denote as shift 0. Next, we double 
the detector size to 2do, and obtain centroid data at the 
two shift positions and do- Analogously, for detectors 
of size to do, we calculate centroid data for to suitable 
shifts. This way, the wave function is estimated at the 
same set of positions for all detector sizes. 

For this analysis, we employ two different methods for 
the calculations for different shifts. First, we generate 
No realizations of the correlated A^-photon state as de- 
scribed in Sec. Ill C 11 In method I, we then use the total 
number of A^o events for each of the m required shifts. 
In method II, we divide the total number of events A^o 
by the number of shifts to, and evaluate each shift with 
No/m events only. In terms of an experimental realiza- 
tion, in method I, the number of measurements required 
increases with detector size, whereas in method II, the 
number of measurements is independent of the detector 
size. 

Due to the statistical nature of the measurement, 
in principle, also a limited number of individual non- 
adjacent detectors can be used, if their positions are 
shifted such that the combinations of all measurements 
cover the entire beam correlation area. This method has 
been adopted in [T3, ■ While the measurement time in- 
creases, the complicated fabrication process for an array 
of detectors can be avoided in proof-of-principle measure- 
ments. 



E. Error estimation 

To evaluate the performance of the centroid method 
for given detection parameters, we first fit the centroid 
data to the original distribution in order to obtain the 
optimum overall scaling factor for the centroid data, and 
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FIG. 3. (Color online) Result of the centroid method for 
NOON states with different N, shown against the centroid 
position in units of the single-photon transverse wavelength 
A. Continuous blue, dashed red, and dotted green curves 
correspond respectively to JV = 2, 3 and 4. The increase in the 
number of fringes clearly shows the resolution enhancement 
with 1/N. 



multiply the centroid data with this factor. This is neces- 
sary, as we evaluate the centroid data only over a limited 
position range, such that the overall normalization of the 
measured centroid data a priori is unknown. We then 
employ the weighted root mean square deviation given 
by 



1 



(33) 



where b is the number of centroid data points, and 
|^(.Xi)| 2 and Zi are the reference value from the orig- 
inal probability distribution and the estimate obtained 
from the centroid method, respectively, at positions JQ. 

Note that the rms deviation in part depends on the 
spatial extent over which the centroid data is compared 
to the original wave function. If a larger range is consid- 
ered, positions are included into the analysis for which 
both the original probability distribution and the cen- 
troid prediction are very low, such that they in essence 
do not increase the sum in the rms deviation, but only b. 
This is particularly important if rms values for different 
wave functions arc compared. 



III. RESULTS 

In the following, we present our results for three dif- 
ferent trial states. We start with NOON states. 



A. NOON States 

1. Resolution enhancement 

Fig.[3]dcpicts the the centroid probability distributions 
for two, three and four photon NOON states, obtained 
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FIG. 4. (Color online) Root-mean-square deviation of the 
recovered two-photon centroid distribution from the original 
wave function position space distribution as a function of the 
detector size. Note the logarithmic scale. 



with a detector size small enough to resolve all features 
of the wave function. It can be seen that there are 27V 
fringes per wavelength for N photons, confirming the ex- 
pected scaling of the obtained resolution with 1/N. We 
thus find that our simulation technique is able to recover 
the predicted resolution enhancement. 



2. Dependence on the detector size 

We thus turn to a detailed analysis of the centroid tech- 
nique. As our first step, Fig. 0] shows the dependence of 
the rms deviation of the recovered two-photon centroid 
distribution from the original position space distribution 
as a function of the detector size. As expected, the de- 
viation is low for small detectors, and initially increases 
with growing detector size. Starting from detector sizes 
of about A/2, the rms deviation saturates. Next to this 
general structure, an increase of the deviation at detector 
size A/4 and 3A/4, as well as sudden drops in the rms for 
A/2 and A can be observed. Finally, towards very low 
detector sizes, the rms deviation starts to increase again. 
In the following, we will explain and interpret all of these 
features. 

First, we analyze the structures at A/4, A/2 and A. 
For this, we choose a detector size, and then calculate 
the rms deviation as a function of a shift of the whole 
detector array along the centroid axis without modifying 
the detector size. The results are shown in Fig. [5] It can 
be seen that for detector sizes other than A/4, A/2 and 
A, the rms deviation is approximately constant over the 
whole range of shifts. For these detector sizes, we con- 
clude that the naive expectation that a smaller detector 
leads to a better recovery of the wave function is correct. 

In contrast, at detector size A/2, the rms deviation 
strongly depends on the shift. Even more surprisingly, 
while there is a slight dependence on the detector shift for 
detector sizes A and A/2, the rms deviation remains very 
low for all shifts despite the large detector size. We will 
show now that these result are artifacts of the calculation 
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FIG. 5. (Color online) Root-mean-square deviation of the 
recovered two-photon centroid distribution from the original 
wave function position space distribution as a function of the 
detector shift. The different curves show results for the detec- 
tor sizes (i) triangles faced right A, (ii) triangles faced left A/2, 
(iii) blue blobs 0.3A, (iv) green squares A/4, and (v) magenta 
rhombuses A/ 10. Note the logarithmic scale. 



procedure, which arise due to the particular structure of 
the wave function to be recovered. 

Fig. [6] explains the situation for detector size A/4. In 
(a), the original position distribution along the centroid 
direction for a two-photon NOON state is shown together 
with the centroid data obtained for detector size A/4 
and zero shift. For this shift value, one of the centroid 
data points coincides with the position distribution max- 
imum at X = 0. It can be seen that half of the centroid 
data points coincide perfectly with the original distribu- 
tion - but the other half strongly deviates. The reason 
is that the position distribution is zero only at single 
points, whereas the large detectors cover a range around 
these zeros with non-zero photo detection probability. 
As a result, the centroid data is non-zero, in contrast 
to the probability distribution. This deviation for half 
of the points explains the large rms value at zero shift 
in Fig. [5j We next consider the same setting, but with 
shift s = — 13A/200, corresponding to a minimum in the 
rms value for detector size A/4 in Fig. [5J Fig.HJb) shows 
the corresponding centroid result obtained with the same 
scaling of the centroid data as in (a) as blue dots. It can 
be seen that the deviation to the probability distribu- 
tion is strong. But if we fit the overall scaling factor of 
the centroid data to the original distribution, then good 
agreement is obtained. This is the reason for the mini- 
mum in Fig. [5J in which the centroid data is fitted to the 
original distribution for each shift position individually. 
The analogous analysis for detector size A is shown in 
Fig. [7] In this case, due to the A-periodicity of the wave 
function, a perfect fit of the centroid data is possible for 
any detector shift. 

However, in an actual measurement, results from all 
shifts would have to be accounted for with the same over- 
all prefactor, as an individual fitting to the unknown dis- 
tribution to be determined is impossible. For this reason, 
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FIG. 6. (Color online) Centroid analysis for detector size A/4. 
The continuous red curve shows the probability distribution 
for a 2-photon NOON state. The blue dots indicate the results 
of a centroid measurement with particular detector settings. 
In (a), the detection array is not shifted with respect to the 
origin of the probability distribution, such that one of the 
centroid data points coincides with the center position X = 0. 
(b) The detection array is shifted to s = — 13A/200, which 
corresponds to a minimum in the rms deviation in Fig. [5] 
The blue dots show the centroid data with the same scaling 
as in (a). The green asterisks show this data fitted to the 
original distribution. 



detectors of size A or A/ 2 cannot recover the wave func- 
tion. To show this, we have taken a set of Nq = 10 6 de- 
tection events, calculated the centroid distributions for a 
large number of shifts, and joined the corresponding dis- 
tributions for the respective shifts without further pro- 
cessing them beforehand. Only the total set obtained 
from all shifts was then fitted in amplitude to the orig- 
inal position distribution. The result was a structure- 
less Gaussian which did not contain the sub-wavelength 
fringe features of the original distribution. Therefore, in 
conclusion, it is not possible to recover the wave function 
with detectors of specific larger sizes, despite the small 
rms deviation for specific detector shifts, which has to 
be interpreted as an artifact of the analysis. Note, how- 
ever, that such measurements with larger detectors could 
potentially be sufficient to analyze the periodicity or the 
symmetry of a given wave function. 

Next, we analyze why the rms deviation increases again 
at very small detector sizes. We found that this is a sta- 
tistical effect. If the detector size is decreased while the 
detection area is kept constant, the number of detectors 
increases. Then, the mean number of events per detec- 
tor decreases with the detector size. At some point, the 
statistical fluctuations due to the decreasing number of 
events become large enough to dominate the rms devia- 
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FIG. 7. (Color online) Centroid analysis for detector size A. 
The figure is analogous to Fig. [6] except for the detector size. 
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FIG. 8. (Color online) Effect of statistical fluctuations on the 
centroid methods for small detector sizes. The rms deviation 
was evaluated (i) black blobs, once for 10 6 events, (ii) red 
squares, twice for two subsets with 5 x I0 5 events, (iii) blue 
rhombuses, 5 times for subsets of 2 x 10 5 events, and (iv) green 
triangles, 10 times for subsets of 10 events. In all cases, the 
same 10 6 events were analyzed, and the respective subsets 
were chosen disjunct. The figure shows the results averaged 
over the respective subsets. Note the logarithmic scale. 



tion. In order to verify this interpretation, we evaluated 
the rms deviation (i) once for 10 6 events, (ii) twice for 
two subsets with 5 x 10 5 events, (iii) 5 times for subsets 
of 2 x 10 5 events, and (iv) 10 times for subsets of 10 5 
events. In all cases, the same 10 6 events were analyzed, 
and the respective subsets were chosen disjunct. After 
this, we averaged the results over the respective subsets, 
such that in every case, the same events were analyzed. 
The result is shown in Fig. [5] It can be seen that the 
rms deviation of all cases approximately agrees towards 
larger detector sizes. But towards smaller detector sizes, 
the smaller the corresponding subsets are, the more the 
deviation increases. Note that also the difference in rms 
deviation between different subsets of equal size is much 
smaller than the difference in deviation between subsets 
of different size. 



3. Comparison of the two methods 

In this section, we compare the two analysis methods 
I and II introduced in Sec. Ill Dl Fig. [3] shows results for a 
2-photon NOON state using the two methods. The two 
curves for method I (blue squares) and II (red dots) agree 
for larger detectors, but deviate for small detectors. In 
particular, the unexpected increase in the rms deviation 
towards low detector sizes observed in Fig. [8] and inter- 
preted as statistical fluctuations due to a low number of 
events per detection bin is absent for method II. 

The reason for this qualitative difference is as follows. 
For the lowest detector size do, no shift is required, and 
the two methods are equivalent. Thus, the same rms de- 
viation is obtained. As can be seen from the increase in 
the rms value for method I, the prediction at this detec- 



tor size is already limited by the low number of events 
per detection bin. For the next higher detector size 2do- 
in method I, the number of events per detection bin is in- 
creased, as the same number of N events are distributed 
over a lower number of bins. In contrast, for method II, 
the number of events per bin is not increased, as with 
increasing detector size, the number of events considered 
for each detector shift is reduced. As a consequence, 
the rms deviation for method II monotonously increases 
with detector size, whereas that for method I can de- 
crease with increasing detector size. In the latter case, 
the reduction of the rms due more events per bin on av- 
erage outweighs the increase in rms due to the increase 
of the detector size. 

This interpretation is further supported by the two 
other data sets in Fig. [SJ These show method II, but 
with half the number of events (green diamonds) or one 
third of the events (black triangles). It can be seen that 
the rms results for the smallest detector size do roughly 
agree with those for method II with half the data points 
at detector size 2 do. Further, the rms for method II with 
one third of the data points agrees to this value at detec- 
tor size 3 do- This shows that for such small detectors, 
the rms deviation is dominated by the low number of 
counts per bin. 

We obtained qualitatively similar results also for N <E 
{3,4}. 



4. Multiphoton NOON states 

We now turn to NOON states with higher photon num- 
ber. Fig. [TU1 shows the rms deviation as a function of the 
detector size for N — 2, 3 and 4 photons. Note that 
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FIG. 9. (Color online) (a) Root-mean-square deviation as a 
function of detector size. Results are shown for a two-photon 
NOON state. The blue squares are obtained using method I 
discussed in the text, whereas the red dots are obtained using 
method II. The green diamonds and the black triangles are 
obtained using method II with one half and one third of the 
total number of events used in (a), respectively. Inset shows 
a magnification at small detector sizes. 
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FIG. 10. (Color online) Root-mean-square deviation versus 
detector size for NOON states using method I. Results are 
compared for number of photons N = 2 (lower blue curve), 
N = 3 (middle red curve), and N — 4 (upper green curve). 
The lines indicate linear fits to the curves over the indicated 
detector size range. 



Fig. [TU] has been generated including all detector shifts 
required for the respective detector sizes, such that the 
spurious maxima or minima in the rms deviation at spe- 
cific detectors sizes found in Fig.@]due to the dependence 
on the shift position shown in Fig. [5] do not appear. All 
data points are generated using the same No = 10 6 events 
with method I introduced in Sec. IIII A 31 As the spatial 
extent over which the probability distribution is signifi- 
cantly larger than zero decreases with increasing photon 
number N, we adjust the position range over which the 
wave function is matched to the centroid measurement 
accordingly. Thus, the 3- and 4-photon cases are eval- 
uated over 2/3 and 2/4 = 1/2 the range of that of the 
N = 2 case, respectively. 

It can be seen that at small detector sizes, the rms 
deviation is low and approximately independent of N. 
In this limit, the detectors are chosen small enough to 
recover all features even of the N = 4 wave function. 
With increasing detector size, the rms deviation starts 
differing from the small-detector limit value first for the 
N = 4, then for the N = 3, and finally for the N = 2 case. 
Increasing the detector size further, the rms deviation has 
an approximately linear dependence of the detector size, 
until it eventually saturates for the large-detector limit. 

In the near-linear region, we estimate slopes a = 
A(rms deviation)/ A(detector size) for N = 2, 3 and 4 as 
0.57, 0.72 and 0.87, as indicated in Fig. 1101 respectively. 
These values are consistent with a scaling with s/N. One 
possible interpretation of this scaling could be as follows. 
If the detection positions of each of the N photons ac- 
quires an uncertainty of order do due to the finite detector 
size, then an estimate of the uncertainty of the centroid 

coordinate X — -hYli=i x i i s gi vcn by SX = d /\/N. 
This estimate is motivated by the fact that the sum of N 
normally distributed independent variables with widths 



8x is again a normal distribution, with width V N8x. Ap- 
proximately, it also holds for other distributions. To- 
gether with the prefactor 1 /N in the centroid variable, a 
scaling of 1/\N in the centroid coordinate uncertainty 
is obtained. Next, we use that an increase of the detector 
size from the low-size limit leads to most significant con- 
tributions to the rms deviation in the regions of highest 
slope of the centroid probability distribution. We denote 
such a centroid coordinate with highest slope as x max , 
and note that the value of the cos 2 () part in the proba- 
bility distribution evaluates to 1/2 at this point. We thus 
can estimate the scaling of the rms deviation via 




cos 2 (27rAa;) - - 



dx . 



(34) 



To leading order in the detector size, this expression 
scales as s/N do, i.e., linearly with the detector size, and 
with slope proportional to s/N, as observed in the nu- 
merical data. 



5. Single- and multi-photon detection 

One motivation for the centroid measurement method 
is the technical difficulty to achieve multi-photon or 
photon- number resolving detection, as required, e.g., 
for sub- wavelength correlated multiphoton measurement 
schemes. If the size of the individual detectors in the 
centroid method is small enough, the number of events 
in which two or more photons hit the same detector is 
negligibly small, such that no photon number resolution 
is required. 

Motivated by this, we have analyzed the percentage 
of two-photon events in our numerical calculations, as a 
function of detector size. We denote the position range 
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FIG. 11. (Color online) Root-mean-square deviation versus 
detector size for a jointly Gaussian state using method I. The 
parameters are same as Fig. {TJ. 
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over which the photon probability distribution is evalu- 
ated by p. In the numerical calculations, p is chosen such 
that it encompasses all parts of the probability distribu- 
tion which are significantly larger than zero. For detector 
size do, we then estimate the number of detection bins p 
as integer closest to p/do- In our numerical simulations, 
we have used p = 7A for N — 2. For N > 2, we used 
accordingly p = 7 A x 2/N. 

We then count those events in which all photons arrive 
in the same bin as a multi-photon detection event. There 
are p different possibilities to realize a A^-photon detec- 
tion event. Overall, in our numerical calculation, there 
are p N different possible events, as we do not make use 
of the symmetrization of the photon states and therefore 
distinguish between events, e.g., in which two photons 
are detected at positions Xi,X2 and X2-,x\. Thus, we 
expect a ratio of multi-photon events given by approxi- 
mately p/p N ■ In terms of the detector size, we find that 
p/p N ~d»-\ 

We then numerically generate random detection 
events, discretize them to model the detection procedure, 
and count all events in which all photons arrive in the 
same detection bin as multi-photon detection events. Fi- 
nally, we calculate the percentage of the such obtained 
multi-photon events out of all detection events. Our nu- 
merical analysis indeed confirms the analytical estimate. 

Note that for the NOON state, the ratio of two-photon 
events reaches about 14% probability for detectors of size 
A. This maximum value may appear high, but it arises 
from the fact that in our calculations the parameters are 
such that the non-negligible support of the probability 
distribution only covers few single-photon wavelengths. 
Then, a detector that is wavelength sized already covers 
a significant part of the total relevant centroid coordi- 
nate range. It is important to note that our calculations 
remain valid despite the two-photon events, since our nu- 
merical approach naturally is capable of photon-number 
resolving "measurements" such multiphoton events do 
not have to be discarded. 
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FIG. 12. (Color online) = A" 1 , (a-c) Total multi photon 
absorption rate and (d-f) peak multi photon absorption rate 
versus the spot size reduction factor for jointly Gaussian state. 
(Top, middle and bottom row) correspond to N (E {2,3,4}, 
respectively. The data points obtained via numerical calcula- 
tion are blue while red is the theoretical prediction. We mul- 
tiplied the Rtot values with 17.8, 156,4100 for N G {2,3,4}, 
respectively to normalize the data such that at r = 1 the 
corresponding value on the y-axis becomes unity too. 



B. Jointly Gaussian States 



We now turn to the analysis of jointly Gaussian states, 
as suggested in [lfl ■ The analysis is performed in a sim- 
ilar way as for the NOON states. First results for the 
dependence of the rms deviation on the detector size in- 
corporating all the shifts at larger detector sizes using 
method I are shown in Fig. (jTTJ) . As for the NOON 
state, the rms deviation increases with detector size, as 
naively expected. The increase of rms deviation towards 
the smallest considered detector sizes is again due to sta- 
tistical effects, as explained in Sec. IIII A 21 Compared to 
the NOON states, the absolute values of the rms deriva- 
tion arc smaller, which is due to the much simpler struc- 
ture of the Gaussian states. As for the NOON states, we 
also analyzed the rms deviation versus different detector 
sizes without taking into account shifts of the detector 
array. Contrary to the case of NOON states, there is 
no qualitative difference between the results with and 
those without shift. The reason is that for the Gaussian 
states, there is no fringe pattern in centroid direction 
which could match the periodicity of the detection array. 
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1. Resolution Enhancement versus Multi photon Absorption 

Next to the accuracy of the wave function recovery, 
also the efficiency is of relevance in any practical imple- 
mentation. In essence, the measurement time is limited. 
This raises the question, how a limited number of mea- 
surements can be used in the most efficient way. In par- 
ticular, in [l8j . the efficiency of centroid detection was 
compared to that of a multiphoton detection scheme. To 
analyze this in our numerical calculation, a practical way 
to distinguish single- from multi-photon detection events 
is required. We thus introduce a distance (Imp = A/400, 
and denote N photons which would hit a detector with 
unlimited position resolution at distance smaller than or 
equal to (Imp a TV-photon event. Photons separated by 
more than (Imp are registered as individual single-photon 
events. 

We start by analyzing the multi-photon absorption 
rate as a function of the spot size reduction factor r, 
which has been defined in [l8[ as the ratio of the classical 
width Wc to the width W of the probability distribution 
of the joint Gaussian defined by Eq. (j2"0")l for different 
choices of the parameters B and /3 such that (fc 2 ) has 
a fixed value. 



0.015 



Explicitly, r = y/NB/^/Jk 2 ). r = 1 de- 
fines the standard quantum limit, whereas the ultimate 
Heisenbcrg limit is given by r = y/N. 

Our numerical scheme works as follows. Wc first fix 
to a positive integer multiple of A -1 . Then, we vary 
B and calculate (3 from Eq. (fill]) . Note that if (fc 2 ) is 
fixed to a large positive value in units of A -1 , then many 
values of r in the quantum regime given by 1 < r < y/N 
cannot be accessed because increasing B results in a /3 
which makes B less than its classical value P/y/N . After- 
wards, we generate random detection events as explained 
m Sec. Ill C II Note that we do not apply discretization to 
mimic detection by finite size detectors. We then deter- 
mine those events which contain photons with distance 
smaller than or equal to (Imp and call these close events. 
The width of the histogram of these close events provides 
us with the width W, which is calculated by fitting the 
histogram data to a Gaussian distribution of the form 
ce~ dx2 . The width W is then given by 2/( sqrt2d). Af- 
terwards, we divide the number of close events by the 
total number of A-photon events to calculate the nor- 
malized total multiphoton absorption rate Rtot- 

Fig. (fT2"k) shows our results for the normalized total 
multiphoton absorption rate function of the spot 

size reduction factor for iV e {2,3,4}. It can be seen 
that the results quantitatively agree with the theoreti- 
cal prediction R tot = ((N - r 2 )/(A - l))(W-i)/2 found 
in [18| . Note that the statistical fluctuation of our nu- 
merical data increases with r towards the limiting value 
y/N. The reason for this is that once the variance of k n 
is fixed, an increase in B is balanced by a decrease in j3. 
At r = Rtot = lj the standard quantum limit is obtained 
which corresponds to a unit multi-photon absorption and 
a classical distribution. As B is increased, the probabil- 
ity distribution becomes more narrow along the centroid 
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FIG. 13. (Color online) Root-mean-square deviation plotted 
versus the detector size for a fixed physical size of the quantum 
mechanical jointly Gaussian centroid probability distribution 
for various N using method I. The blue, red and green curves 
correspond respectively to TV = 2, 3 and 4. Fhe parameters 
are explained in the text. 



direction, but expands along the orthogonal directions of 
the relative coordinates. This implies that the distance 
between any two random coordinates can be larger now. 
This is why with an increase in B, the number of close 
events decreases. For this calculation, we extracted 1 
million pairs of random events. At B = 0.01 A -1 , out of 
the 1 million events, some 79, 201 are close events while 
at B = 0.999A -1 , this number reduces to 3316. 

Using our r, we can calculate normalized peak multi 



photon absorption rate defined in [18[ as r((N — r )/(N 



l))^ -1 )/ 2 . Fig. (fT2"t O shows the corresponding result. 
The blue data points show the numerical results while 
red curve plots ry/2 — r 2 . Again, the two results agree 
at most of the range of spot size reduction factor r but 
r = y/N is not approached by the numerical data. 

Note that we obtained qualitatively similar results also 
for N = 3 and N = 4 photon states, shown by Fig. (|12b'l. 
(e) and (c), (f), respectively. 



2. Resolution Enhancement versus Multi photon 
Absorption for a Fixed Feature Size 

So far, the analysis has been limited to cases in which 
the feature size decreases with increasing number of de- 
tected photons N. Now, we compare the rms deviation as 
a function of detector size for different N, but with phys- 
ical feature size kept constant. For this, we notice that 
the probability distribution in centroid direction scales as 
cxp(— 2N 2 B 2 X 2 ), with X the centroid coordinate. For 
each A, we choose B in such a way that the resulting spa- 
tial extent of the distribution becomes exp(— 8(X/X) 2 ), 
which leads to B = 2/ (NX). We further choose (3 = 1/A 
for N = 2,3. For N = 4, we instead choose /3 = 4/5 /A, 
such that in all cases the parameters fulfill B > (3/ y/N 
and thus the probability distribution is non-classical. 

The results presented in Fig. ([TBI show that as N is 
increased, the value of rms deviation for a fixed detector 
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FIG. 14. (Color online) Root-mean-square deviation against 
detector size for a correlated coherent cat state with a — i. 
Results have been obtained using method I. 
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FIG. 15. (Color online) Root-mean-square deviation versus 
detector size for a correlated coherent cat state with a = i. In 
contrast to Fig. 1141 results for a single fixed detector position 
are shown. 



size is decreased for most of the range of detector sizes. 
This is to be expected because if the physical size of the 
distribution is fixed, absorbing more photons should im- 
prove the result of the measurement. As explained be- 
fore, the results at very small detector sizes are prone 
to numerical fluctuations if a fixed number of events is 
distributed over more and more detector bins. 

Qualitatively similar results have also been obtained 
using method II. Again we found that method II leads 
to slightly higher rms deviations in particular at lower 
detector sizes, because of the lower number of events in- 
cluded in the calculation. 



C. Correlated Coherent Cat States 
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FIG. 16. (Color online) (a) Probability distribution for a cor- 
related coherent cat state with a — iy2 in position space, 
(b) Root-mean-square deviation against detector size using 
method I. 



Finally, we turn to coherent cat states. As for the other 
states, we start by discussing the rms deviation against 
the detector size. Results are shown in Fig. [H] Note 
that this figure was obtained including all different shift 
positions for the larger detectors, using method I such 
that the positions at which centroid data is obtained are 
the same for all detector sizes. The results is as expected, 
with a certain range of small detector sizes over which 
the recovery of the wave function is good, followed by 
a continuous increase of the rms with the detector size, 
until the rms saturates towards large detector sizes. 

In contrast, the rms deviation for single fixed detection 
arrays of different sizes show a rather different behavior 
compared to the other considered states, as shown in 
Fig. 1151 Here, an initial increase in rms deviation with 
the detector area is seen, as expected, but for detector 
sizes of about 0.17A, the rms deviation acquires a max- 
imum and afterwards oscillates with the detector size. 
These features again can be explained via the shape of 
the probability distribution of the original wave function. 
Results for other shifts of the detection array give similar 
results for detector sizes up to about A/10, but afterwards 
deviate significantly. 

As can be seen from Eq. (f2~T)) . cat states have the ad- 
vantage that the fringe pattern in centroid direction su- 



perimposed onto the overall Gaussian wave function en- 
velope has a periodicity which can be controlled by the 
magnitude of a. Since |ct| 2 is the mean number of pho- 
tons in the corresponding cat state, this change in the 
fringe pattern goes along with a change in the number of 
photons in the light field. But nevertheless, some of the 
realizations of the cat state will consist of two photons 
independent of \a\, such that we can continue to evalu- 
ate only those detection events. In this way, the feature 
size can be controlled without changing the number of 
detected photons. An example for a correlated coherent 
cat state with higher mean photon number (a = iy2) is 
shown in Fig. (fTS)) . As compared to Fig. [ljf) with lower 
mean photon number, the number of fringes in the posi- 
tion space probability distribution increased as expected. 
Note that a similar shrinking was observed in Fig. [3] for 
NOON states with higher number of detected photons. 
Here, in contrast, the number of detected photons is kept 
constant at 2. 

Fig. 1171 shows the rms deviation as a function of |a|. 
With increasing \a\, the feature size of the wave function 
to be recovered becomes smaller, such that for a fixed 
detector size, increasing rms deviation with increasing 
| a | is expected. One can notice that the rms deviation 
initially remains approximately independent of the mag- 
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FIG. 17. (Color online) Root-mean-square deviation for cor- 
related cat states as a function of \a\. In all cases, two-photon 
detection events are considered. No = 10° events were con- 
sidered with detector size A/100. 



number of different values for ranging from to ir/2 
and found that the rms deviation as a function of the 
phase has some residual fluctuations which are com- 
patible with a constant value independent of the choice 
of 0. It is interesting to note that the position space prob- 
ability density as a function of the two position variables 
x\ and X2 is rather different for = tt/2 than for = 0. 
For = 0, it consists of two near-circular structures. 
But for = 7r/2, the probability density is strongly elon- 
gated along the direction perpendicular to the centroid 
axis, and narrowly-spaced fringes along the centroid axis. 
Thus we conclude that the centroid method is suitable to 
recover non-periodic patterns along the centroid axis, in- 
dependent of the precise structure of the position space 
probability density along the directions perpendicular to 
the centroid axis. 
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FIG. 18. (Color online) Position space probability density 
along the centroid coordinate for \a\ = 1 and (i) = 0, (ii) 
= tt/8, (hi) 4> = 3708, (iv) 4> = tt/2. 

nitude of \a\, and then increases with growing |a|. We 
can thus conclude that the parameters chosen such that 
for |a| below approximately 2.5 the detector size is small 
enough to recover the wave function as good as possible, 
limited by the statistical uncertainty due to the finite 
number of detection events. Higher \a\ lead to feature 
sizes which cannot be recovered fully. 

Finally, we exploit that the coherent cat states allow 
one to continuously tune the shape of the centroid po- 
sition space probability density from a periodic fringe 
pattern to less periodic structures. This can be achieved, 
e.g., by choosing different phases 0, as shown in Fig. IT51 
For = 0, an approximate double-peak structure is 
achieved. = 7r/2 leads to the multi- fringe pattern 
studied so far, with relative heights of the peaks given 
by the overall envelope of the wave function. But for 
intermediate values, the shape of the wave function can 
be changed, such that, for example, the central peak is 
suppressed (0 = tt/8. This way, the shape of the prob- 
ability density can be modified from the periodic fringe 
pattern observed, e.g., for NOON states. This is signif- 
icant not least due to the dependence of the rms devia- 
tion on the structure of the wave function we observed 
in Sec. IIII A 21 Wc performed the centroid analysis for a 



IV. SUMMARY 

In summary, wc have analyzed the optical centroid 
method for superresolution lithography using numerical 
experiments. In the respective limits, our results fully 
agree with the predictions from the previous approximate 
analytical analysis. But our numerical treatment allows 
us to gain additional insight, also beyond the validity 
range of the analytical predictions. We exploited this 
mainly by studying the different obscrvablcs in detail as 
a function of the detector size. A first analysis revealed 
unexpected features in the root-mcan-squarc deviation of 
the recovered wave function compared to the true wave 
function. In the limit of low detector size, the rms does 
not continuously decrease as expected, but eventually in- 
creases again. This could be traced back to fluctuations 
in the number of events per detector bin which grows 
with decreasing detector size if the number of events is 
kept constant. In addition, a fixed detection array can 
lead to spurious maxima or minima in the rms deviation 
at specific detector sizes. In particular, we found very low 
value of rms deviation for large detectors of certain sizes. 
These features are found to originate from the structure 
of the wave function, but in general do not allow to re- 
cover the wave function with high resolution using large 
detectors. Such effects, however, could be used to detect 
certain properties of wave functions such as symmetries, 
or to efficiently distinguish between different wave func- 
tions using an array of large detectors. We then aug- 
mented this analysis with results for NOON states with 
higher photon number, and with results for jointly Gaus- 
sian states with fixed feature size, but different photon 
number. Interestingly, for the NOON states, we found 
a regime of intermediate detector sizes in which the rms 
deviation depends approximately linearly on the detector 
size. One plausible interpretation of this effect is in terms 
of the error in estimating the centroid coordinate arising 
from the non-zero detector size, evaluted at the part of 
the wave function with highest slope which contributes 
most to the rms deviation. In the final part, we showed 
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that the OCM can also be applied to cat states made 
out of superpositions of coherent states. The coherent 
cat states have the interesting property that their spatial 
feature size can be continuously tuned via the magnitude 
of the coherent state parameter a. Also the structure of 
the cat state position space probability density along the 
ccntroid axis can be modified via the phase of a. This 
way, for example, shapes other than the periodic fringe 
pattern observed, e.g., for NOON states, can be probed. 
We found that the centroid method works equally well 
for different structures of the cat state. Since the num- 
ber of photons in the cat state input pulse is not fixed, 
we focused on detection events with a certain number of 



photons. In this sense, our results in this section also 
relate to post-selection methods. 
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